Entanglement in Dilute Flux Line Liquids 
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We construct, within the framework of classical statistical mechanics, a mean field theory of 
dilute flux line liquids which goes beyond linear hydrodynamics. Within our approach, we find that 
interactions between vortices produce a massive term in the Hamiltonian of the internal modes of 
the flux lines which confines their transverse fluctuations. This suggests that the flux line liquid, 
at least in the low density limit considered in this paper, might very well be in a weakly entangled 
state, where the average width (u 2 ) 5 of flux lines can be much larger than the average distance 
between the lines but does not diverge with the sample thickness L. Consequences on the physics 
of flux line liquids are briefly discussed. 
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I. INTRODUCTION 

The discovery, by Bednorz and Mulleifl in 1986, of 
the new family of cuprate-based high temperature su- 
perconductors (HTSC), has had a profound impact on 
research in the field of superconductivity. In partic- 
ular, due to their anisotropy, short coherence lengths 
and high critical temperatures, it quickly became obvi- 
ous that thermal fluctuations will play a key role in the 
physics of interacting flux lines in these materials. Early 
experimentsoij -suggested, and it was later unambigu- 
ously confirmed0~tL3, that a phase transition takes place 
where the three dimensional AbrikosovEil flux line lattice 
(FLL) lattice melts into a flux line liquid, and the loca- 
tion of the melting temperature T m as a function of the 
applied magnetic field H in the (H , T) phase diagram of 
the flux line system became the subject of intense experi- 
m&fttal and theoreticaL-H-- investigations (see Blatter et 
alxl for a recent review) . 

One of the early approaches aimed at describing the 
melting transition and the physics of melted flux line liq- 
uids was the theorvpOf-Nelson and Seungt3, which is based 
on the observatior£3'tHl that the partition function of a 
system of interacting flux lines in 2 + 1 dimensions can 
be mapped onto the imaginary time partition function 
of a system of interacting quantum mechanical bosons in 
two dimensions. This formal mapping was used to ob- 
tain a complete description of flux line liquids by using 
well known techniques from boson physics. The phys- 
ical picture of the flux line liquid which emerges from 
this description suggests that there is considerable wan- 
dering of flux lines as they traverse the sample and that 
the flux line liquid is heavily entangled over much of the 
liquid phase. Such a picture has been observed in nu- 
merical simulations, and it is now generally believed that 
the mean square relative displacement ([r^(z) — ri(0)] 2 ) 
of the transverse position r^(z) of a given flux line in a 
flux line liquid has the same form as the corresponding 
quantity for an isolated flux line, namely : 



the only effect of the interactions between lines in a dense 
liquid being to reduceEfl the "diffusion" constant D from 
its, bare value D for a single flux line, which is given 
byBB : 



D 



k B T 
£l 



where ks is Boltzmann's constant, T is the temperature 
and i\ is the tilt modulus of a single flux line. Similarly, 
the mean square projected area occupied by thn-LLnes is 
found, within the boson picture, to be of ordeiMtJ 
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(1.2) 



where u(z) = r(z) — (r(z)) and L is the sample thickness 
in the direction of the flux lines (which we take parallel 
to the z axis). This last result is to be contrasted with 
the mean square projected area of apgiven flux line in the 
crystalline phase, which is given byEj : 



(u 2 ) = 



ksTa 



(1.3) 



(|r l (z)-r,(0)| 2 )=2D|z| 



(1.1) 



where a is the average spacing between flux lines in the 
vortex lattice, e = A/A c is the ratio of the London pen- 
etration depths in the (ah) plane and along the direc- 
tion of the c axis respectively (we remind the reader that 
most HTSC have uniaxial symmetry, and we denote by 
c the principal axis of symmetry of the crystal), apd 
eo = (0o/47rA) 2 , with (j)Q = hc/ 2e t he flux quantumES. 
Comparing the results (1.2) and (1.3), we see that, while 
the mean square projected area of flux lines (u 2 ) in the 
crystalline phase is finite, in the liquid phase it diverges 
with the sample thickness L. n n 

Despite a good_.d£aL.of both theoreticat__E_l and ex- 
perimepiaLwork__r__j'__j and rather extensive numerical 
studiesL-JE-l, several questions regarding the physics of 
flux line liquids remain open. One of the most im- 
portant questions which has been the--sjftbject of in- 
tense debate during the past few yeara__li--l, is related 
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to whether a disentangled flux line liquid can survive 
above melting. Although early numerical simulations 
have |-Hrovided support for the proposal by Feigel'man 
ct a/.Ej that superconducting coherence could survive in 



the liqui 
authorr' 28 



more extensive simulations by several 



found that longitudinal coherence van- 
ishes simultaneously with melting, and that the vortex 
liquid is characterized by very short entanglement lengths 
in the direction of the flux lines. However, analyses of 
experiments on high quality untwinned-single crystals pt 
YBa 2 Cu 3 7 _ 5 (YBCO) by Righi et alB and by MooreEll 
have suggested that longitudinal correlations in the liquid 
phase may be of surprizingly larger micron scale. 

At the fundamental level, it is very important to un- 
derstand, from a purely conceptual point of view, under 
what circumstances a disentangled state can emerge from 
a theoretical description, and what kind of theories are 
best suited to yield such a disentangled state. We first 
note that, for a disentangled state to appear, it is nec- 
essary that flux line conformation variab les vAz) have 
correlations which do not obey equations (1.1)-(1.2). In 
particular, the average projected area of a given flux line 
(u 2 ) has to be independent of the sample thickness L. In 
the crystalline phase, transverse fluctuations of a given 
flux line are confined thanks to the repulsive action of 
the other, neighboring lines in the FLL. In a similar way, 
physical intuition suggests that the interactions between 
flux lines should strongly reducecj and might even sup- 
press line wandering if the repulsion energy between flux 
lines is high enough compared to kgT. In such a dis- 
entangled state with large repulsion forces between the 
lines, a given flux line will experience the repulsive poten- 
tial of its surrounding neighbors in much the same way 
as in a lattice. As a result, the internal fluctuations of 
that line will be confined in about the same way they are 
confined in regular Abrikosov FLL. In particular, there 
should be a regime where the mean square internal fluctu- 
ations (u 2 ) does not differ much, at least locally, from the 
corresponding quantity in a FLL at the same density. For 
us to be able to describe a putative disentangled phase 
of a flux line liquid, that is, a phase where fluctuations of 
flux lines are confined, we need to concentrate on the in- 
dividual conformation variables of the flux lines, and try 
to evaluate the effect of interactions on these variables in 
the most accurate way possible. 

In fact, most theoretical studies to date have concen- 
trated on hydrodynamic quantities, such as, e.g. the den- 
sity, and have used approximations pertinent to linear 
hydrodynamics in order to study the physics of flux line 
liquids. Hydrodynamic observables being coarse-grained 
quantities, such approaches are intrinsically unable to 
give an accurate picture of the correlations of flux line 
conformation variables, such as Ui(z). EvexL-if we go be- 
yond linear hydrodynamics, such as in refo, because of 
the use of coarse-grained variables, we expect not to be 
able to describe the confinement of individual flux line 
fluctuations that would take place in a hypothetical, dis- 
entangled flux line liquid. To the most, hydrodynamic 



approaches can only yield a renormalization of the tilt 
modulus of flux lines in a line liquid, leaving the quali- 
tative (analytic) form of the correlations ([Ll])-([L^) un- 
changed. 

In view of the above remarks, it seems to us that it is 
highly desirable to construct a theory which keeps track 
in a better way of conformation variables of the flux lines, 
and which would therefore be able to yield more accu- 
rate information about the correlations of these variables. 
Shuch a theory would undoubtedly help us gain a better 
understanding of the properties of flux line liquids. 

In this paper, which is largely motivated by the con- 
clusions of referenceEy, we would like to lay out the path 
for the construction of such a theory. We show, in partic- 
ular, that a classical theory of flux liquids with a better 
handling of the internal modes of flux lines can be formu- 
lated which is simple enough to allow for ease of math- 
ematical treatment, and which is transparent enough to 
allow for approximations to be made with a good level 
of confidence. Within this approach, and under certain 
assumptions, we find that a confinement of the internal 
fluctuations of the lines can indeed take place, at least 
in the limit of a dilute flux line liquid. More specifically, 
we find that, due to the interactions between flux lines, 
the internal modes may actually become massive, leading 
to a finite mean square width (it 2 ) 2 of flux lines, and to 
a qualitative change in the behaviour of the correlation 



function (f.l) 



This paper is organized as follows. In section we 
consider a liquid of interacting flux lines in a weakly 
anisotropic HTSC. Performing a perturbative expansion 
of the interaction energy in terms of the internal fluc- 
tuations of the flux lines, we show that the interactions 
with neighboring lines induce a mass term in the Hamil- 
tonian of the internal modes. As a result, we find that 
the internal fluctuations of the flux lines are strongly sup- 
pressed compared with the free flux line result (|L^) . We 
self-consistently determine the range of validity of our 
perturbative expansion, which is found to be quite large 
for moderately anisotropic superconductors. We also cal- 
culate the structure factor of the flux line liquid and our 



result is compared to previous derivations. Section III 
contains a discussion of our findings, which we think can 
be reconciled with the results of numerical simulations, 
along with our conclusions. 



II. MEAN-FIELD THEORY OF (PUTATIVE) 
DISENTANGLED FLUX LINE LIQUIDS 



In this section, we shall be considering the statistical 
mechanics of an assembly of flux lines in a sample of 
thickness L of a uniaxial HTSC, with both the princi- 
pal axis of anisotropy c and the external magnetic field 
aligned with the z axis. Such a system can be described 
by the HamiltoniantL3 
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(2.1) 



where £\ ~ £ 2 £ l n (a/£) denotes the tilt modulus per 
unit length of the flux lines (here £ and a are the co- 
herence length and the average distance between flux 
lines in the (ab) plane respectively), and where V(r) = 
2eoKo(r / A) is the interaction potential between flux line 
elements— at, equal height, with Kq a modified Bessel 
functionE3o. (Henceforth, we shall neglect the logarith- 
mic factor ln(a/£), of order unity, in the definition of the 
tilt modulus, writing i% w £ 2 £o-) The Hamiltonian above 
may be viewed as a simplified version of the most general 
Hamiltonian of arbitrarily curved and tilted flux lines in. 
the London regime of anisotropic superconductorsEiT lj 
(with nonlocal interactions in the z direction), and can 
be shown to be a good approximatipa to the latter in 
the limit of nearly straight flux linesE3. The tilt modu- 
lus per unit length E\ on the other hand can be viewed 
as the effective tilt modulus for long-wavelength q z dis- 
tortions which accounts for the most relevant fluctuations 
of the flux lines near melting, with transverse wavelength 
<Z_l ~ Ibz (siBZ being the transverse wavevector at the 
Brillouin zone boundary), as can be seen by replacing 
q_i_ ~ qsz = (47r_B/(/)o) 1 / 2 , taking the limit q z — > in 
the general (nonlocal) expression of the tilt modulus per 
unit volume C4 4 (q) = B 2 /4tt(1 + (\/e) 2 q\ + A 2 g 2 ), and 
multiplying by the elementary area per vortex (cf>o/B). 
Below, it will prove useful to write the following Fourier 
decomposition of T{(z) into Rouse modest : 



v i (z)= J2 r Mj qnZ 



(2.2) 



where q n = 2irn/L, and where the coefficients r, (g n ) are 
given by : 



r i(<Zn) = 7" / dz r,(z)e lq " 



(2.3) 



as can be verified by using the orthogonality relation 

f-L 



dz e l «" z (e l « mZ )* =LS„ 



(2.4) 



It will also be convenient to write i"j(z) as the sum 

Vi(z) = r 0i + VLi(z) (2.5) 

where 

roi = ri(q n =0) = — I dz r^z) 
L Jo 

is the center of mass (CM) position, while Ui(z) is the dis- 
placement of the flux line at height z with respect to the 



center of mass position. In terms of the Fourier modes 
{ri(g„)}, the displacement vector Uj(z) can be written in 
the form (here c.c. denotes complex conjugation) : 



u i( z ) = ^2{ri{q n )e lqnZ + c.c.} 



(2.6) 



Using the above decomposition of the flux line posi- 
tions {ri(z)} into CM and internal modes, equation (glj), 
the interaction term between flux line elements at i"j(z) 
and Yj(z) can be written as 

V(v l {z) - r 3 {z)) = V(r Ql - r 0j + u,(z) - u,(z)) (2.7) 

Since the interaction potential V(r) is a smooth function 
which varies very slowly on length scales much smaller 
than A, we see that if 

ImOO-uj-coKA (2.8) 

then we can write, to a very good approximation : 

V(r i (z)-r j (z))~V(v oi -r oj ) + 

+ [Ui(z) - Uj(z)} • W(r 0i - r 0j ) + 

+ ±([ui(z) - u;(z)] ■ V)V(r w - r oi ) (2.9) 

If we integrate this last equation over z, taking into ac- 
count the fact that J Q dz Uj(z) = 0, we obtain 



dz V(vi(z) - Tj(z)) ~ LV(r 0i - r oj ) + 



dz [u ia (z)Uip(z) +U ja (z)Ujp(z) + 

- u ia {z)u jl3 {z) - u ja (z)u i0 (z)] d a dpV(r Qi - r 0j ) 

In the Boson language, this last equation is reminiscent of 
the expansion of-the action of a quantu m pa rticle around 
a classical pathO. The Hamiltonian (2.1) now can be 
written in the form : 



H = Hn 



Hi 



where 



1/0 = ^1%-^) 



(2.10) 



(2.11) 



is the Hamiltonian of a system of perfectly straight, in- 
teracting lines of length L, and 



N 



dz 



1 _ /du 2 \2 1 (i ) 

£ 1 ( — ) + T, V a p u ia{Z) U if 3{z) 



V dz 



i=l 

^EE/ dz n { $u ia {z)u jf) {z) (2.12 
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is the Hamiltonian of the internal modes of the flux lines. 



In equation ( 2.12 ), we defined the following quantities 



Vaf! = -d a dpV(r Qt -r j) 



(2.13) 
(2.14) 



As can be seen from equation ( p. 12 ), the interac- 
tions between flux lines have generated a mass term 
^fj£gUi a (z)uip(z) for the internal modes {ui(z)} as well 

as additional, two-body terms, \lJ>^a Uia(z)ujp(z) , which 

couple the internal modes of different vortices (i =/= j). 

fii (a) 

As they stand, the coefficients p af} and pr a p of equations 



(2.13) and ( |2. 14 ) depend on the specific configuration of 
all the center of mass coordinates {r oi } considered. In 
order to obtain the value of these coefficients relev ant to 
a liq uid of flux lines, we need to average equations ( 2.13| ) 
and (2.14) over all possible configurations of the CM co- 
ordinates compatible with a liquid structure. In order 
to carry out such an average, let us write the partition 
function of the flux line liquid (/3 = 1/fceT is the inverse 
temperature) in the form 



Z = Tr Tr u (e 



-0H\ 



(2.15) 



where we denote by TVo and Tr u the trace over the cen- 
ter of mass and internal modes, respectively. More ex- 
plicitely : 



Tr 



dr i • ■ ■ dr ( 



" 



(2.16) 



Tr„ = / II dri ^ / II dr N(Qn) (2.17) 

n>l "* n>l 

where dri(q n ) stands for dY^ re (q n ) dr i>im (q n ); r hre (q n ) 
and Yi^ m {q n ) being the real and imaginary parts of 
Ti(q n ), respectively. Noting that H = Ha + H\ in ( |2.15 ), 
and taking the trace over the CM mode only, we obtain 



Tr (e-^)=Z (e-^>o 



(2.18) 



where Zq = Tro(e~@ Ha ), and where the average (■ ■ -)o 
is taken with respect to the probability distribution 
exp(— (3Hq)/Zq, i.e. for an arbitrary function /({i"oi}) 
of the CM coordinates {r i, • • • , r Ar}, 

</({r 0l })) = ±-Jdr 1 ...dr N f({r 0l }) e"^ (2.19) 



We n ow perform a cumulant expansion in equation 
( 2.18 ), with the following result to leading order in Hi, 



Tr (e-^) ~ Z e 



-/3(Hi) 



(2.20) 



and hence we see that the total partition function Z can 
be written in the form 



ZqZi 



(2.21) 



where Z\ = Tr u (e P( Hl )°) is the effective partition func- 
tion of the internal modes. The quantity 



H, 



eff 



(Hi)o 



may be thought of as the effective Hamiltonian of the in- 
ternal modes, averaged over all possible CM liquid con- 
figurations. In order to find this effective Hamiltonian, 
all we need to do is to calculate the averages (/i„n)o an d 



(fJ-ag)a- Applying the definition of the average ( 2. IE ) to 



(i) 



the "mass" coefficients p a/3 



we obtain 



(A 



i r N N 

T / IT dr ° k E d adpV{r Qi - r 0j )e~ 



0H O 



N 

E 



\- \ dr 0l dr Q:j d a dpV(r 0l - r oj ) 
J 



- I Y{ dr ok e~ 



I3H 



(2.22) 



The coordinates {tqi} under the integral sign bein g jus t 
dummy variables, we see that the rhs of equation ( 2.22| ) 
consists of a sum of (N — 1) identical terms which can be 
rewritten in the form 



(mS)o 



(N 



- 1) f 

; / dr 01 dr 02 d a dpV{r 01 - r 02 ) x 

o J 



-< J dr 03 . . . dr 0N e 

i J dr 01 dr Q2 d a dpV(r 01 - r 02 ) 
N(N - 1) 



dr, 



■ dr 0N e 



(i Ho 



(2.23) 



We recognize on the second line of equation ( |2.2 
two-particle density pW(ri, r 2 ) of liquid state thcor 



P 



(2) (roi,r 02 )^ N{N _ 1] /,/r„ ,...,/,•, 



Z 



on e" 



I3H (2.24) 



This last quantity is related to the pair distribution func- 
tion of the CM mode <7o(roi,ro 2 ) by 

p (2) (r i,r 02 ) = p 2 #o(roi,r 2) 

where p = B/cjjQ is the average density of flux lines in 
the system. At equilibrium, we expect the flux line liq- 
uid to be translationally invariant, which implies that 
go( r oij r 02) will depend only on the difference (r i — r 02 ), 
i.e. 3o( r oij r 02) = 5o( r oi — 1*02) - Under these conditions, 
we obtain from equation ( |2.23 ) 



2 /> 

(^)o = Jf dr oidr 02 ffo(roi - r O2 )d a d V(r Q1 - r 02 ) 
= p [dr 9o (r)d a dpV(r) (2.25) 
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Introducing the pair correlation function ho (r) = go (r) — 
1, the rhs of the last equation becomes 



g HrL 



(0\ 

a/3/0 



p dr h (r)d a di3V{r) 



(2.26) 



where we used the fact that J dr d a dpV(r) = 0. If we 
further use the fact that both ho and V are rotationally 
invariant functions, h {r) = ho(r) and V(r) = V(r) , we 
obtain (here d± = 2 is the number of transverse dimen- 
sions) 



(MS)o 



P 5 a ,p 



P = 



where = d 2 



f [dr h(r)V\V{r) 



(2.27) 
(2.28) 



d 2 is th e Laplacian in the transverse 



directions. Equation ( 2.28 ) is the main result of this pa- 
per. It gives the effective mass of the internal modes of 
flux lines, averaged over all possible configurations of the 
centers of mass of the v ortices in the flux line liquid. 

From equation ( 2.28 ) above, we now want to derive 
an explicit expression for the averaged mass p. In or- 
der to be able to do so, we need to choose a suitable 
analytic form for the pair correlation function ho(r). It 
is well known from liquid state theoryEPJO that the pair 
distribution function g(r) of an ordinary classical liquid 
has peaks of decreasing amplitude at r rs a, 2a, 3a, . . .; 
a = 1/ ./jo being the average spacing between particle in 
the liquid phase; and that it has the limiting behaviors 
g(r) <C 1 for r < a, and g(r) — > 1 when r — > oo. A com- 
mon (mean field) approximation for the pair distribution 
function used in the literature consists in neglecting the 
fine structure of g(r) for r > a, and approximating 



g(r) ~ 9(r - a) 



(2.29) 



where 6 is the Heaviside unit step function (6(x) = 1 if 
x > and zero otherwise) . This approximation, which is 
obviously relevant to liquids at not very high densities, 
retains the important feature of g(r) that it is very small 
for r < a, which expresses the fact that particles in a liq- 
uid of density p have a small probability of being within 
a distance a = 1/sfp from each other. Here, we expect 
such a behaviour of the pair distribution function go(r) 
in a dilute flux liquid to be valid if the repulsive interac- 
tions between flux lines are strong enough. To quantify 
how strong the interactions must be, we introduce the 
following typical energy scale 



E C (T) = eo l c (T) 



(2.30) 



where l c (T) = max(s, £ C (T)), s being the distance be- 
tween superconducting planes in our layered material. 
E c represents the typical energy barrier for the crossing 
of two flux lines nearly parallel to c, and if E C {T) ksT 
we expect interactions between lines to strongly reduce 
cutting and crossing of lines, and as a consequence, to 




Figur e 1. Pair distribution function go(r) of equation 
(2.31) for q = 3/2. Also shown is the mean-field ap- 
proximation go(r) — 0(r — a). 

lead to the kind of confinement of the internal modes 
we find above. In this regime, equation ( 2.29] ) should 
be a good approximation to the pair distribution func- 
tion go (r) in a dilute flux liquid. Although the condition 
E c 3> ksT may prove to be too stringent for actual flux 
lines in a HTSC, we here consider this limit as an extreme 
case for which a disentangled state might be realized. 

For the purpose of mathe matic al tractability, we shall 
here use the approximation ( 2.29| ) in the form (figure 1) 



g (r) = 1 - exp(-w 2 /a 2 ) 
which gives for the pair correlation function 
ho{r) = — exp(— ar 2 /a 2 ) 



(2.31) 



(2.32) 



The numerical coefficient a inside the exponential is to 
be chosen so that \ho(r = a)| < 1. Usi ng thi s form of the 
pair correlation function into equation (2.28), we obtain : 



/' : 



d±a 



(pa 2 ) 



k 2 V(k)e- k2a2 / 2a (2.33) 

[2-kY 



Inserting the Fourier transform of the interaction poten- 
tial 



V(k) 



47T£r 



k 2 + A- 2 

we obtain, after a few manipulations (Appendix A) 



(2.34) 



2tt 

p w — peo 

d±_ 



(2.35) 



where we assumed that a < X. It is interesting to note 
that the final result is, to a very good approximati on, in - 
dependent of the particular choice of a in equation (2.32), 
as long as this last condition (a < A) is satisfied. It is also 
interesting to note that a nontrivial structure in go{r) 
is necessary to obtain a finite value for the "mass" p : 
inserting in equation ( 2.25| ) the trivial (hydrodynamic) 
approximation 



9o(r) = 1 



< r < oo 



would lead to a vanishing result for p and to the super- 
fluid kind of behaviour described by equations ( 1 . 1 )- ( 1 .2 ) 
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We now turn our attention to the off-diagonal terms We thus obtain for Gij(q n ) the following result 



(ij) 

H a p which couple the internal modes of vortices i and j. 
We have : 



(M, 



«/3 /0 



iV 

fc=i 

--7T ( dv m dr 02 d a dpV(r 01 - r 02 ) x 
< f dr a3 . . . dr 0N e~ 0H ° 



Comparin g the rhs of the last equation with the rhs of 
equation (|2.23|), we see that 



1 



(^oia ) - (JV- 1) 



(2.36) 



which, in view of the fact that the mass \i = {fjy') is fi- 
nite (i.e. independent of N), implies that the off diagonal 
terms vanish in the thermodynamic limit N — > oo. 

We are now left with the following averaged Hamilto- 
nian for the internal degrees of freedom of our flux line 
liquid : 



N 



dz 



1 _ /rfuA 2 1 2 



H eff = E 

-W^J)Y,Y,[ dz^u^-u^z) (2.37) 



N 



where we used the results (2.27)-(2.36) for the averaged 



diagonal (/i^g) and off-diagonal terms, in conjunc- 



tion with equation ( 2.12| ). This last expression of H e ff 
can be rewritten in Fourier space in the form 

N 1 

H °ff = E E 2 i^'h^rM ■ rj (- 9n ) (2.38) 

i=l n^O 

with the inverse propagator 
(G~%(g„) = [Lei^+i/i + 



iV-1 



N — 1 



From this last equation, the propagator Gij(q n ) defined 
through 



(ui(q n )uj(q m )} = k B TG(q n )S n 



(2.39) 



can be easily obtained using an identity for inverting 
N x N matrices of the form 



namely, 



(A 1 ) ij = aSij 



A - — —A 



a(a + bN) 



(2.40) 



(2.41) 



ij[(ln) - T(q n ) ° ij + T(q n ) [T(q n ) - LpN/(N 1)] 



where we denote by T(q n ) the quantity 

L/i 



r(q n )= [Le x ql + Ln- 



N - 1 



(2.42) 



(2.43) 



In th e ther modynamic limit N — > oo, we see from equa- 
tion ( 2.42 ) that Gij(q n ) reduces to the diagonal form 



Gij(q n ) = G(q n )Sij, with 
G(q n ) 



Le x In + Lfx 



(2.44) 



Hence, the total projected area of a given flux line in the 
liquid environment is given by 

(u 2 (z)) =d ± k B TY,G(<ln) 



2d ± k B T 
Lex 



y L . 

^ Ql + (A*/ei) 



(2.45) 



Using the result ( p.35| ) in equation ( 2.45 ), along with the 
fact that 



oo 

El 7T , . I 

—z = — coth(7ra) 

n 2 + a 2 2a 



n=X 



2 a 



2a 2 



for a — > oo 



we obtain (we take c/j_ = 2) 

/ 2\ k B T 



2ee JiTp 



(2.46) 



(2.47) 



The important thing to note about this result is that, 
unlike the free flux line, the mean square projected area 
(it ) is not proportional to and in fact does not depend 
at all on the sample thickness L. More significantly, if we 
take p ~ 1/a 2 , where a is the average intervortex distance 
in the liquid phase, the above equation gives us : 



k B T a 
2^Tree 



(2.48) 



which is exactly t he r esult one obtains in the crystalline 
phase, equation ( jl.3| ), and is actually what we expect 
on purely physical grounds for a flux line trapped in the 
cage potential formed by its surrounding neighbors, as- 
suming that the barriers for cutting and crossing are high 
enough. That we are able to obtain such a result for (u 2 ) 
gives us confidence in our approach and in the various 
approximations that were made in the course of our cal- 
culation. 
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Having obtained (u 2 ), we are now in a position to as- 
sess rtiie range of validity of our analysis. Using the fact 
thatB s (K/A) = 1.964 x 10 8 /[A(A)] 2 , and the values 
e = 1/7, A(0) = 14001 relevant to YBCO, we find 



where p(r, z) is the local density operator at height z 



N 



(2.54) 



u 2 ) ^ 0.27 x 10~ 2 f A(0)^ 

2 ~ 



a 



1 - T/T c 



(^) 



(2.49) 



where we approximated J~p ~ 1/a, and where we used 
the mean field dependence of the London penetration 
depth on temperature A(T) = A(0)/^1 - T/T c . The 
width of the flux lines in a dilute liquid (a < A) is hence 
very small compared to the average inter-vortex distance 
a. We need to go to much higher densities (i.e. to a 
higher ratio A/a) and to temperatures very close to T c in 
order to have (u 2 ) « a 2 . More quantitatively, in order for 
our perturbative analysis based on the taylor expansion 
( |2.9D to be correct, the condition ( pT§| ) has to be satis- 
fied. This condition can be rewritten in the following, 
averaged form : 



According to the theory of referenceEl, the partial Fourier 
transform S*(q, z) should behave as 



S(q,z)~S(q,0)e 



-e(q)\z\ 



(2.55) 



where e(q) is the Bogoliubov spectrum of the correspond- 
ing bosons. The above form of S*(q, z) suggests that the 
densities at two different heights z and z' become more 
decorrelated as \z — z'\ grows. Here, due to the fact that 
the internal modes are massive, we expect the density 
to remain correlated on the whole longitudinal length of 
the sample. Indeed, in Appendix B we show that the 
structure factor S(q, z) within our model is given by : 



([u l (z)-u J (z)] 2 )<a 2 



for all ijtj (2.50) 5(q, z) = p 2 5 o(q) e" 



But, since (iij(z) ■ Uj(z)) = (11,(2)) • (uj(z)) 



for 



i ^= j in our approach (we remind the reader that the 
off-diagonal coefficients py^) were found to vanish in the 
thermodynamic limit), the condition above can be rewrit- 
ten in the form 



+ pexp - q 



2y/neeo 



1 



•i*h 



(2.56) 



(2.51) 



which, by equation (2.49), is sure to be satisfied on a 
rather large region of the phase diagram (H, T) of the 
flux line system. 

Let us now find the difference correlation function 
([vLi(z) - Ui(0)] 2 ), which is given by : 

, 2 4d±k B T ^ 1 - cosq n (z - z') 
([u(z) - u(z )] > = 2^ 



Le~i 



n=l 



Transforming the sum into an integral, and using the re- 
sult 



dq 1 " C ° S( f } = f [1 
q 2 + a 2 2a 1 



-a|x|] 



we obtain, after a few manipulations : 



({u(z)-u(z')} 2 } 



d±ksTa 



\z-z'\\ 



a > 



(2.52) 



which shows that the structure factor consists of two 
pieces : a CM piece, proportional to p 2 , which corre- 
sponds to a liquid of more or less straight flux lines and 
includes, through ga(q L ), nontrivial correlations between 
the CM positions of flux lines, and a second piece linear in 
the density p, which describes the internal fluctuations of 
the individual lines and which incorporates the effect of 
interactions between vortices through the confining mass 
/1. Like the structure factor of the boson analogy, our 
5(q, z) does decrease as a function of \z\, with the im- 
portant difference, however, that 5(q, z) here goes to a 
finite limit (much like in a crystal} as \ z\ grows very large, 
in contrast to the result of ref.E3 where 5(q, z) — * as 
\z\ — > 00, indicating that flux lines remain correlated over 
rather long length scales in a disentangled flux line liquid. 

The most important implication of the theory devel- 
opped in this section is related to the issue of entangle- 
ment of flux lines. Assuming that the typical energy scale 
for the interactions between flux lines is much higher than 
ksT and neglecting thermally nucleated vortex loops, we 
have derived the statistical mechanics of a dilute line liq- 
uid using a simple Taylor expansion of the Hamiltonian 
in terms of the internal modes of the flux lines. We have 
calculated the mean square projected area of a given flux 



We thus find that the mean square displacement ([u(z) - linc ( u )_ and we have found that is of thc same OTder 
u(z')] 2 ) goes to a finite limit as \z — z'\ — * 00. This 



was expected, since in our approach, thc internal modes 
are massive and have bounded fluctuations which do not 
grow with the sample thickness. 

The last quantity we shall be interested in is the 
density-density correlation function 



S{v-v' ] z-z') = {p{v,z)p{v',z')) 



(2.53) 



of magnitude as what is usually obtained in a flux line 
lattice. The physical picture of a flux line in a liquid envi- 
ronment that emerges from our analysis, under the above 
assumptions, is that of a roughly straight object whose 
internal modes are contained in a tube of finite radius 
(u 2 ) < a, no matter how thick the sample might be. At 
this stage, we can already assert that we have achieved 
the goal stated in the Introduction, and that we now have 
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at our disposal a detailed and yet simple theory which is 
able to describe a completely disentangled flux line liquid. 
As we mentioned earlier, the assumption E c ksT is, 
however, too strong for real flux liquids on the melting 
line H m (T) (as can be easily checked using material pa- 
rameters relevant to most HTSC), and our theory needs 
to be adjusted to take this fact into account. In the fol- 
lowing section, we discuss and compare our findings with 
previous work, and we try to construct a consistent pic- 
ture of what we expect the physics of flux line liquids to 
be based on our approach. 



III. DISCUSSION AND CONCLUSIONS 

We shall start this section by briefly reviewing and dis- 
cussing the results of numerical simulations of flux line 
liquids. A number of these simulations have used the 
mapping of the Ginzburg-Landau free energy of a super- 
conductor onto the Hamiltonian of the uniformly frus- 
trated three dimensional (3D) XY model : 



Hxy 



cos(#j — 9j — Ai 



(3.1) 



where, in a discrete lattice, 0, is the phase of the 
complex Ginzburg-Landau order parameter ip{ r y z ) — 
\iP(r,z)\e l9( - r ^ at site i, A tJ = (27r/0 o ) // A • dl is the 
line integral of the magnetic vector potential along a 
path linking site i to site j, and the sum is over nearest- 
neighbor sites. The Hamiltonian (3.1) results from mak- 
ing the London approximation ipQ = ^(x, z)| = Cst. 
in the discretized Ginzburg-Landau free energy of an 
anisotropic superconductor in the regime A/a — > oo, so 
that the magnetic induction B = Bz inside the super- 
conductor can be taken as a constant. For a layered su- 
perconductor with the principal axis c along the z axis, 
the couplings Jjj are given by Jjj = J± when sites i and 
j arc located on the same superconducting plane, with 
J± = ^s/(16tt 3 A 2 ), and J ifj = J z = £0^ 2 /(16vr 3 sA 2 ) 
when sites i and j are located at ±sz from each other, 
where s is the spacing between superconducting layers. 

Within the model ( |3.1[ ), evidence for an entangled state 
is deduced from various dynamical quantities, such as the 
helicity modulus T(q) which gives the linear response of 
the supercurrent j to a transverse perturbation in the 
vector potential A ext of the externally applied magnetic 
field : 



(3.2) 



The helicity modulus measures_the. response of the sys- 
tem to an imposed phase twistlZjiZj, and is therefore re- 
lated to the correlation function (here R = (r, z)) 



*(R) = (V(R)^*(0)expi 



2l7T 



»o Jo 



R 



A- dl 



(3.3) 



which measures the coherence of the phase degrees of 
freedom in the system. Also measured in numerical sim- 
ulations of model (3.1) is the structure factor 5(r, z) — 
(p(r, z)p(0, 0)), where the density operator p(r,z) = 

12iLi ^( r — T i( z )) i s given in terms of the phase variables 
0(r, z) by 

p(x,z) = ^[Vx V0(x,js)] 'Z 

Information about entanglement and longitudinal corre- 
lations is then deduced from the behaviour of the struc- 
ture factor 5(0,2;) = (p(r, z)p(r, 0)) at zero transverse 
separations. 

Another model which has been used in numerical stud- 
ies of flux liquids is the so-called lattiee._Lpndon model, 
which is defined by the HamiltonianOE^EJ 



* = 2*VlE E 



G^TLi - R^MR^n^R,-) (3.4) 



where n M (Rt) is the integer vorticity through plaquette 
p at site Ri = (r^z^) of a cubic mesh of points, and 
J_l is the same coupling as the one defined above in the 
context of the uniformly frustrated 3D XY model. The 
lattice London interactions G M have Fourier transforms : 



G,(q) 



Ql 



Ql 



(Ql 



2 /A 2 



(Q 2 + e 2 S 2 /A 2 )/(g 2 + . s 2 /A 2 ) 



Ql + Q 



{Ql 



VA 2 ) 



(3.5) 



where — 2sin(g M s/2) and Q 2 = Tl^Q 2 ^- Unlike the 
3D XY model, which has Coulomb like interactions be- 
tween vortex elementscH and which is only valid in the 
regime A 3> a of dense systems, the lattice London model 
takes screening into account and can therefore be used for 
dense as well as for dilute flux liquids. 

The last kind of simulation which has been per- 
formed. is. the Monte Carlo simulation ofpNordborg and 
BlatterBO who use the boson mapping)! 3 ] to simulate 
a system of interacting flux lines. In this approach, 
the authors adapt an algorithm which has been used in 
the pastLj to study superfluidity in quantum-mechanical 
Bose systems, involving large scale permutations of vor- 
tex trajectories subject to the "periodic" boundary con- 
ditions rj(0) = Tj(L), i.e. every line ends either on itself 
or on another line. Below we comment on this Monte 
Carlo method and the associated algorithm, but, before 
doing so, we want to consi der t he two other simulation 
methods, models (3.1) and (|3.4|). 

We first observe that a direct comparison between our 
results and sim ulations of the uniformly frustrated 3D 
XY model ( |3.l| ) is made difficult by the fact that, while 
this model is relevant to the case A 3> a, we here con- 
sider the opposite limit a < A of a dilute flux line liquid. 
In addition, and as we already mentioned, most of the 
simulations above were concerned with the measurement 
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of the helicity modulus, which is essentially a measure 
of phase correlations in the system. It should be noted, 
however, that phase correlations can be very different 
from density correlations..— To illustrate this statement, 
we refer the reader to ref.113, where it was found within 
a simple elastic ap pro ach that the order parameter for 
phase correlations (^^) decays exponentially in the direc- 
tion of the superconducting planes in a three dimensional 
flux line latticed, although the density order parameter 
Sk(t,z) — (pic( r , ^)Pk(O) 0)) is known to be finite in 
such a lattice. This example shows clearly that there is 
no "one to one" correspondence between phase and den- 
sity correlations, and that the vanishing of the helicity 
modulus in a flux line liquid for example does not nec- 
essarily mean exponential decay of density correlations 
along the z axis or correlations be t ween flux line posi- 
tions of the type of equations This is even 
more so if we allow for the formation of vortex loonsL. 
Thermally excited vortex loop excitations were foundE£l 
to destroy phase coherence and to lead to a vanishing 
of the helicity modulus—even in an otherwise perfectly 
ordered flux line latticed 

Finally, it should be realized that the simulations of 
models (|3 . 1[) and (3.4) are highly simplified representa- 



tions of considerably more complex physics. For example, 
simulations are usually carried out at a given density of 
flux quanta / = (Bt; 2 /</>o) = (B /2irH C 2) (we remind the 
reader that the upper critical field for H||c is given by 
H C 2 = 0o/27r£ 2 ), with values of / ranging from 1/12 to 
1/30 most frequently used. Values of / in this range 
correspond to values of the induction B in the range 
~ \H C 2 — \Hc2. Apart from the fact that the use of 



models (3.1) and (3.4), which are based on the assump- 
tion of a constant amplitude of the order parameter, be- 
comes questionable in this range of fields, we also note 
that at such high fields vortices in adjacent superconduct- 
ing planes4n the liquid phase might very well be already 
decoupled!!! (note also that the temperature dependence 
of the coherence length £ in / = (B^ 2 /<f>o) and of the Lon- 
don penetration depth A in J± and J z is generally left out 
as T is varied). Moreover, in most of these simulations, 
which are carried out on a cubic mesh, the mesh constant 
in the z direction is associated with the distance s be- 
tween superconducting planes in a layered material, and 
no mention is made whatsoever of the coherence length 
£ c along the c axis. In particular, the question regarding 
whether the same results, e.g. for density correlations in 
the liquid phase, would follow in the regimes £, C (T) < s 
and £, C (T) > s, has largely remained untouched. The 
fact that most of these simulations consider implicitely 
the quasi- two-dimensional case&J £ C (T) < s does not al- 
low us to draw any conclusions on the interesting regime 
£ C (T) 3> s where the average tilt angle of flux linesE 2 ] 
{{dvi/dz) 2 ) ~ (kBT/e 2 eo£, c (T)) can be very small and 
where we expect superconducting coherence, whether it 
be for the phase of the superconducting order parame- 
ter or for the density, to survive on much longer length 
scales. 



Several remarks are due regarding_the Monte Carlo 
simulation o£-Mordborg and BlatteiE3'Ej who use the bo- 
son mappinglid to simulate a system of interacting flux 
lines. First, the use of the "periodic" boundary condi- 
tions Yi{L) = rj(0) imposes additional constraints on the 
system which are not present in a real flux line liquid. Ac- 
tually, in order to "capture the effects of Bose statistics" , 
the flux lines in this simulation were made to switch their 
endpoints by hand, and this might very well introduce an 
artificial entanglement in the system. More specifically, 
the switching of endpoints is achieved by cutting out suf- 
ficiently long segments of a number of lines and trying dif- 
ferent ways of connecting the loose-ends while satisfying 
the periodic boundary conditional. While this proce- 
dure and the corresponding algorithm are appropriate for 
the imaginary time paths of quantum mechanical bosons, 
we find it more realistic that the flux line system should 
be allowed to equilibrate with free (as opposed to peri- 
odic) boundary conditions and, more importantly, with- 
out explicit switching of endpoints. Within the context 
of quantum mechanical bosons, superfluidity is brought 
about in numerical simulations by precisely these kind of 
manipulations. In the context of vortices, we expect these 
manipulations to lead to an overestimation of the effect 
of entanglement in flux liquids, and the "superfluid" be- 
haviour found in referenceo might therefore be a direct 
consequence of the specific bosonic algorithm used in this 
simulation. 

In view of all the above remarks, and the results of 
section |], it seems to us that there is still room, at least 
in the dilute limit considered in this paper and neglect- 
ing vortex loops, for a new phase of the flux line liquid 
which has not been considered in the past, and which we 
might describe as weakly entangled. By "weakly entan- 
gled" flux liquid, we mean a liquid phase in which {u 2 ) 
can be very large, but does not actually diverge with 
the sample thickness L. For this to happen, a small 
value of the confining mass fx, much smaller than the 
one found in equation ( 2.35 ), is needed. In the following 
paragraph, we try to understand how such small values 
of fi can emerge from our model, and how a weakly en- 
tangled phase can be reconciled with what is observed in 
numerical simulations. 

The key quantity in the derivation of the confining 
mass of section || is the pair distribution function </o( r ) 
of the CM mode. In equation (2.31), we used an approx- 



imation for go(r) which is relevant to a situation where 
the average distance between the centers of mass of the 
flux lines is a — 1/y/p, and where the CM positions are 
strongly anti-correlated for r < a, i.e. if the CM of the 
ith flux line is at location roi, then the centers of mass 
of neighboring vortices have a very small probability to 
be within a distance a from r oi . While this is a perfectly 
legitimate way of thinking for actual flux line elements 
and for the regime E C (T) ksT considered in section 
||, since the CM is only a mathematical construct on 
one hand, and since the condition E C {T) ksT is usu- 
ally not satisfied near melting on another (making situa- 
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Figure 2. Pair distribution function go{r) of equation (3.6) 
for a = 3/2 and n = 0.1 (upper curve). For compari- 
son, we also show the pair correlation function of equation 
(2.31), (lower curve). 

tions where CM positions of different flux lines are very 
close to each other rather frequent) , the CM distribution 
function does not have to be very small (go{r) ^ 1) for 
r < a. However, because of the repulsion between flux 
lines, we cannot completely neglect correlations between 
the center of mass positions and simply take go{r) = 1 
for < r < oo , in which case we would r ecov e r a f ully en- 
tangled line liquid with the correlations (1.1)- (1.2). Since 
there is obviously no easy way to derive an analytic ex- 
pression for the pair distribution function go Mi w ^ here 
are bound to speculate about its possible shape using our 
physical intuition and infer the resulting physics. One 
possible shape for go(r) which allows for close encounters 
between the centers of mass of different flux lines is the 
one shown in figure 2, and corresponds to the approxi- 
mation 



.9o (r 



(r) ~ 1 — r\ exp ( — ar 2 /a 2 ) 



(3.6) 



where rj is a numerical constant such that < r\ < 1 (in 
figure 2, we use rj = 0.1). Such a form of go(r) would 
yield a confining mass which is n tim es smaller than the 
one found previously, equation ( 2.35 ), and would there- 
fore lead to values of the average projected area (u 2 ) 



ksT a 
2?7 % /7r£eo 



(3.7) 



which arc much larger than the square of the average dis- 
tance a between flux lines. On the other hand, equation 
(3.6) would lead to a relative displacement 



'M2\ 



<[u(z)-u(z')] 



\z-z'\\ 



(3.8) 



For r\ <C 1 there will be a substantial range of separations 

(3.9) 



z-z <^ = 



ea 



where the rhs of equation (3.8) can very well be approx- 
imated by 



([u(z)-u(z')f)^N- Z 'l 



(3.10) 



(note that factors of rj have cancelled each other in 
this last equation) and for which density correlations 
(p(q, z)p{— q, 0)) will decrease exponentially as a func- 
tion of \z\, leading to the structure factor 



5(q, z) ~ pe 



(3.11) 



Note, in particular that for moderately anisotropic super- 
conductors, and if r/ <§C 1, the length £ z can be quite large, 
maybe even larger than the relatively small thicknesses 
used in numerical simulations, giving the semblance that 
t he e ntan gled state is characterized by the correlations 
( 1 . 1 )-( 1 .2 ) with an entanglement correlation length 



2 2 

e e a 
d± ksT 



On the other hand, for values of n close to unity (rj < 1), 
we recover a flux line liquid where correlations extend 
over much longer distances, which-, corresponds to the 
situation analyzed by Righi et al£Q. We thus see that 
our model can describe, depending on what <?o( r ) m re- 
ality is, very different situations corresponding to nearly 
disentangled (n < 1), weakly entangled (n <C 1) or fully 
entangled (77 = 0, <?o( r ) = 1) h ux hne liquids. 

To summarize, in this paper, using a perturbative ex- 
pansion of the Hamiltonian of interacting flux lines in a 
vortex liquid similar to the expansion of_thc action of a 
quantum particle around a classical patbO, we have con- 
structed a mean field theory of the flux line liquid which, 
in the author's opinion, has a better handling of the in- 
ternal fluctuations of flux lines, and which, through the 
use of a nontrivial pair distribution function go(r) for the 
positions of the centers of mass of flux lines, goes beyond 
Gaussian hydrodynamics, making contact with the stan- 
dard liquid state theory of classical fluids. Within our 
approach, we find that a weakly entangled phase might 
exist, where the average width of flux lines (u 2 ) 2 can be 
much larger than the average intervortex distance a, but 
does not diverge with the sample thickness L. By vary- 
ing the shape of go (r) , we are able to describe situations 
with very short and rather large entanglement correlation 
lengths. Since the slightest deviation of the CM pair dis- 
tribution function from unity at small distances (r < a) 
will necessarily lead to a finite (nonzero) confining mass 
fx, a careful measurement of go(r), for example in numer- 
ical simulations, would provide an unequivocal way of 
verifying the predictions of this paper, and to ascertain 
whet h er a fully entangled liquid phase with correlations 
( 1.1 )-( 1.2 ) or the weakly entangled phase proposed here 
is the correct ground state of the flux line liquid. 
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IV. APPENDIX A : DETAILS OF THE 
CALCULATION OF THE MASS fi FOR THE 
INTERNAL MODES OF VORTICES IN A FLUX 
LINE LIQUID 



7to 2 , . 2a\ 2 
(P e o) ■ — T 



d±aX 2 
2tt 



dj_ 



(4.6) 



In this appendix, we show details of th e calculation which is the result quoted in the text, equation ( 2.35 ). 
of the integral on the rhs of equation (2.33) which gives 
t he v alue of the effective mass /i. Using the expression 
(2.34) of the interaction potential in Fourier space U(k), 
and performing the integration over the polar angle, we 
obtain 



2tt 



(pa 2 ) 



i 2 

k 2 + \~ 2 



(4.1) 



V. APPENDIX B : DETAILS OF THE 
CALCULATION OF THE STRUCTURE FACTOR 
OF A FLUX LINE LIQUID 

In this appendix, we show some details of the calcula- 
tion of the structure factor 



Using the change of variables x — X 2 k 2 , this last expres- 
sion can be cast in the form 



S(r,z;r',z') = (p(r,z)p(r' ) z')) 



(5.1) 



dx 



l + x 



d±aX 2 

with k = a 2 /2a A 2 . Now, using the resulJ 

1 



dx 



l + x 



e K Ei( K ) 



(4.2) 



(4.3) 



where Ei(x) = dte t /t is the exponential integral, 
and using the Taylor expansion for small arguments!^ 



Ei (a;) = —7 — In x — 



(-l) n x r - 



n(nl) 

n=l y ' 

(here 7 = 0.577. . . is Euler's constant), we obtain 



(4.4) 



both for a free flux line liquid of nonintcracting flux lines, 
and for the weakly entangled flux line liquid described in 
the text. In general, we expect flux liquids at equilibrium 
to be translationally invariant, so that the structure fac- 
tor will depend only on the relative coordinates (r — r') 
and (z — z'), i.e. S(r, z; r', z') = S(r — r', z — z'). As 
a consequence, we have for the Fourier transform of the 
density-density correlation function : 

(p(q, dz)M> «£)) = (2^) rf <5(q + q')S(q z + q' z ) S(q, q z ) 

(5.2) 

In the following, we shall be concerned with the quantity 



d±a\ 2 



(pe )- -(l + Ke K (7 + lnK) + o(/t 2 )) (4.5) 



S{q,q z ) = —4-; (p(q, «z)p(-q, -Qz)) 



(5.3) 



where we used equation (5.2) above and the fact that 
Using the fact that k = a 2 /2a\ 2 < 1 in the regime a < X (27r) rf 5(q = 0)S(q z = 0) = LLr^ in the limit L, L±_ — ► 00. 
(assuming that a > 1), we finally obtain We have : 

(p(q,<fe)/3(-q,-<fe)) = I dvdz fdr'dz' (p(r,z)p(r',z')) e - < *('-*V < «*<*-* , > 



N N 



e -iq-(r-r') e -iq z {z-z') 



= J2J2 1 drdz I dr ' dz ' ( 6 ( r - r *( z )) 5 ( r ' - r j( z '))) 

i=l j=l ■> ■> 

= dz J dz ' (^ ,q ' (rorro,) )o(^" 1 ' lu ' ( ^' 1,(z ' ,) )ie~" l(z ^' ) + 

i=l jjii 

+ £/„/^(e-*<".«-»*')>) 1 e-*<,-.0 

1=1 

N 

= ^2^2 dz dz ' (e _,q '' r,rr ° i V~^°'' i ''°''°' z '~' Il °' 2 ''''" i '' i ' z '~" iJ ' z '^ 1 q-^Az-z') _|_ 

z=l jjii 
N 

■!■ ^ d,Z ( fz'e~59° < ?/3([ tl i,= ( 2 )- tl i^( z ')][" l ,^(^)-«i,/3(2')]>l e -"lz(z-z') J-g ^ 



where (• • -)o and (• • -)i denote averages over the CM and internal modes with statistical weights exp(— (3Hq) and 
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exp(—(3H e f f) respectively (Ho and H e ff are the Hamiltonians given in equations (2.11) and (2.37) of the text). Both 
for the liquid of noninteracting lines, and for the weakly entangled liquid that we consider here, the internal degrees 
of freedom u,(z) and Uj(z) belonging to two different lines i ^ j are decoupled. We therefore can write, for i ^ j : 



([ui, a (z) - Uj }0l (z!)}[ui,p(z) - u jt p(z')])i = (ui. a (z)ui^{z) + u^ a (z')uj^(z'))i 

- 2<5a/ V> 



(5.5) 



where, in going from the first to the second line, we used the fact that (ui^ a u i p)i = S a ^(uf a ). Equation ( |5.4| ) becomes 

JV r r 

(p(q,««)p(-q,-&)) =EE / dz / dz ' (e-^- r ^) e-^ q2{u2) e^' 2 ') + 

i=l jyti J 

N 

+ S~] / dz / dz'e~^ 



\q a q[i{['u. i , a (z)~u iia (z l )}[u iil3 (z)-Ui 4i (z')]) 1 e -iq z (z-z') 



(5.6) 



Using the fact that f dzf dz' e lqz<yZ z ' = L 2 S q ^ , and noticing that 

N 



-iqr ; p -Jq ■roj\ n — 



)o = (2^) !i ^(q + q') p 2 30 (q) 



(5.7) 



i=l jjii 

which gives us here (with q' = — q) 

N 



J2J2{c-^e^) = (2n)^5(q = 0) p 2 . 90 (q) = i£Vffo(q) 



(5.8) 



we finally obtain : 



(p(^Qz)p(-^-qz)) = L%„ L{± p 2 ff0 (q) e~^ q2{u2) 



J dz J dz 



> c - 2J^1 2 ([ u i(z)-u,(z')] 2 ) i e -iq z (z-z') 



where we used the fact that ([u iya (z) — Ui ia (z')][ui t p(z) — Ui,p(z')\)\ — ([Uj(z) — Uj(^')] 2 )i. From this point on, 
things will differ depending on whether we consider a liquid of free flux lines or a weakly entangled flux liquid. We 
shall treat both cases separately, starting with a liquid of free flux lines. 



A. Liquid of free flux lines 

In a liquid of free flux lines (or in a fully entangled "superfluid" state, for that matter), the average width of flux 
lines (it 2 ) 2 diverges with the sample thickness L. We can therefore drop the first term on the rhs of equation ( pjj ), 
which contains the factor e d ±- q {w ' . If we in addition use the fact that 

([ Uj ;(z) - U,(z')] 2 )l = \Z - Z'\ (5.10) 

£l 

we obtain : 

(p(q,8.)p(-q,-9.)) = £ / dz / dz'e"^^ lz ~ z ' 1 e"^"*') 

i=l J 

dZe~* q l z l e - tq * z (5.11) 

-oo 

Performing the Z integration, and using the fact that the density p — N/L"^ , we finally obtain 
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q 2 + (k B Tq 2 /2 £l ) 

From the structure factor 5(q, g 2 ) = (p(q, q z )p(— q, — q z ))/LL d ^- , we can deduce the density-density correlation func- 
tion 5(r, z) = (p(r, z)/5(0,0)) in real space. We have : 

d d± q /" dq 



S(t,z) 



The integrations being quite straightforward for the free flux liquid considered here, we only quote the following 
intermediate result : 

S(q )Z )=p exp(-^^ \z\) (5.14) 

along with the final result 

which shows that the density-density correlation function (p(r, z)p(r, z')) at the same transverse location r behaves 
like 

s <°-"-'(wrF^r) (5 - 16) 

and thus we see that, even for the fully entangled noninteracting flux line liquid, density-density correlations 
(p(r, z)p(r, z')) at the same transverse location r decrease only algebraically as a function of the height separa- 
tion \z — z'\. It is therefore very surprizing that S(0, z) was found to decrease exponentially, like ~ gxp(—\z\/£ v ), in 
the presence of interactions in some numerical simulations of the uniformly frustrated 3D XY modelEZI. 

B. Weakly entangled flux liquid 

In a weakly entangled flux line liquid, (it 2 ) does not diverge with the sample thickness L, and therefore the first 



term on the rhs of equation (5.9) has to be kept. Using equation (|3.8|), 



([uf» - u(z')] 2 ) * %=|^(1 ~ e"^ lz - z ' 1 ) 



we obtain 



03(q, fe)p(-q, -9,)) - L\ z , p 2 5o (q) e"^^ + NL J dZ exp ( - q 2 ^gg- (l - e" 1 ^ ^) ) e"^ 
= LLi^LS^opigoiq) e"*^ + p j dZ exp ( - 9 2 ^gg-(l - e"^ "*") ) e"<**} 
from which we see that 5(q, g 2 ) = -^^(^(q, g z )/5(-q, -&)) is given by 

5(q, &) = L^, p 2 30 (q) e^^ + p | exp ( - q 2 ^^(l - e"^ ) (5.17) 

where we used the fact that (u 2 ) w k B T a/2r]^/Trseo) , equation ( |3.7[ ). In view of the fact that the integration on 
the rhs of this last equation cannot be performed exactly, it is more convenient to use the partial Fourier transform : 

5(q, z) = p 2 5o(q) e-« 2 m + pexp ( - q 2 ^g- (l - e"^ H) ) (5.18) 

which, with 77 = 1, is the result (2.56) of section [El]. On the other hand, for 77 <C 1, th e "Debye- Waller" factor 



exp ( — q 2 2 ^^£ £o ) can be very small. Neglecting the first term on the rhs of equation (5.18), we obtain 



s<q..W«p(-«" 5 X;(i-^ !P '-) 



which, for |z| < ^ 2 = eajr\\pK leads to equation (3.11) of the text 
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